Designing Isotropic Interactions for Self- Assembly of Complex Lattices 
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We present a direct method for solving the inverse problem of designing isotropic potentials 
that cause self-assembly into target lattices. Each potential is constructed by matching its energy 
spectrum to the reciprocal representation of the lattice to guarantee that the desired structure is 
a ground state. We use the method to self-assemble complex lattices not previously achieved with 
isotropic potentials, such as a snub square tiling and the kagome lattice. The latter is especially 
interesting because it provides the crucial geometric frustration in several proposed spin liquids. 
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Self-assembly [ij, the process through which simple en- 
tities spontaneously form complex structures, is a phe- 
nomenon of fundamental interest and a compelling tool 
for future nanoscale fabrication. As technology is pushed 
to ever smaller length scales, direct fabrication becomes 
increasingly challenging and material design may turn to- 
wards tailoring interactions which causes constituents to 
order into the desired patterns without directed external 
influence. This vision has been the subject of intense ex- 
perimental work during the last decades but there is still 
a relative lack of theoretical results beyond pure simula- 
tion based work. 

What theoretical work exists can be understood as 
trying to solve an "inverse" statistical mechanics prob- 
lem, finding interactions that causes self-assembly into 
a target structure or configurations with desired prop- 
erties. Various schemes for approaching this problem 
have been proposed [H, Q and several isotropic poten- 
tials self-assembling into low-coordinated lattices have 
been reported, e.g. square, honeycomb, simple cubic, and 
diamond lattices [4]. The focus is on isotropic interac- 
tions, partly due to theoretical tract ability but mainly 
motivated by the recent years' remarkable improvement 
in understanding and design capabilities of such inter- 
actions in especially colloidal and nano systems [HI, 
By combining and tuning forces including, but not lim- 
ited to, van der Waals, electrostatic, magnetic, and steric 
forces as well as forces of purely entropic origin we can to- 
day construct complex effective interactions. With such 
a profusion of possibilities, the forward statical mechan- 
ics' prescription of sweeping the parameter space to give 
a complete picture of the possible structures becomes 
computationally infeasible, and the inverse formulation 
increasingly appealing. 

In this Letter we present a new method for designing 
potentials for targeted self-assembly based on an energy 
spectrum approach. The method is direct in contrast 
to iterative or optimization approaches more commonly 
used to solve inverse problems. We show how interactions 
that cause particles to organize into the kagome lattice 
can be obtained. The kagome lattice is a much sought af- 



ter topology due to its highly frustrated geometry 0, IH , 
possibly resulting in spin liquid ground states in several 
materials [ol, [lo| as well as in the corresponding Heisen- 



berg model . The self-assembly of nanoscale 12, 13| 
and colloidal [1^ kagome lattices through directed inter- 
actions has caused considerable interest but whether this 
target structure is achievable via isotropic potentials have 
until now been an open question (cf. the negative result 
in [15]). To demonstrate the generality of our method, 
we also apply it to the honeycomb lattice and a snub 
square tiling [l6| . 

We consider self-assembly governed by energy mini- 
mization where an interaction is designed with certain 
features so that the target structure is a stable low en- 
ergy state, which should typically be reached when the 
system self-assembles from random initial configurations. 
Our method differs from previous ones in that we use 
a reciprocal representation of the lattice to construct a 
theoretical potential in a way that guarantees the target 
structure to be the unique ground state. The potential 
is then systematically simplified into a more realizable 
form. In general a lattice structure is defined by the 
small |k| region of the Fourier spectrum while the large 
|k| region only contributes with higher resolution of the 
distribution defining the particles. The kagome lattice 
for instance is clearly resolved (Fig. [T]) when expressed 
as the sum of the Fourier modes of the (weighted) re- 
ciprocal lattice vectors up to length 2|kp|, where kp is 
a primitive reciprocal lattice vector. This makes the re- 
ciprocal perspective a natural approach for targeted self- 
assembly of lattice structures, as opposed to real-space 
methods focusing on the pair correlation function [4]. 

In two dimensions the potential energy for a config- 
uration with particle density p{r) can be expressed in 
reciprocal space as 

E = y"dridr2p(ri)F(|ri -r2|)p(r2) 

= Jdk \m\'V{k), (1) 

where \p{k)\'^dk = /j^^^ dk |p(k)p, p(k) is the standard 
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Fourier transform of p(r), and V{k) is the radial Fourier 
(Hankel) transform of V{r) [I^, 

V{k) J rdrV{r)Jo{rk). (2) 

If ]/(ri, r2) = V(|ri — r2|) is interpreted as a linear oper- 
ator, then Eq ([T]) defines a diagonalization of V(ri,r2). 
We therefore refer to V{k) as the energy spectrum of the 
potential [18i] . 

We use the power spectrum of the target lattice \p{k)\ 
and Eq. ([T]) to construct a potential that guarantees that 
the target is a ground state in the following way. Con- 
sider a potential V(r) whose energy spectrum is chosen 
positive, V{k) > 0, and has zeros only at those points 
coinciding with the reciprocal lattice points of the 
Bravais lattice describing the target lattice, y(|G^|) = 
Vi. This design guarantees that the target lattice is a 
ground state with E = 0. However, any lattice with re- 
ciprocal points C G^, i.e. with arbitrary basis (as well 
as reciprocal lattices with opposite chirality), is also a 
ground state. We break this degeneracy in favor of the 
target lattice by introducing small perturbations of the 
energy spectrum at its zeros and use the fact that there 
are amplitude differences in the basis dependent p{Gi). 

To understand how perturbations affect the energy lev- 
els there are two characteristics of p(k) we need to con- 
sider. For a closed system with n particles p(0) = n so 
V{0) is irrelevant to the ground state configuration. More 
importantly, the limited mass provides an upper limit for 
|p(k)| = Ill-Li^**^ '^* I < n, with equality for certain k if 

coincides with the resonance points of the modulation. 
It is this property we use to break the degeneracy. 

We now apply this design scheme to the kagome lat- 
tice. Consider an energy spectrum with zeros corre- 
sponding to the hexagonal Bravais lattice with prim- 
itive vectors {(2,0),(l,\/3)}, which encompasses the 
kagome lattice (basis {(0, 0), (1, 0), (1/2, ^3/2)}). The 




FIG. 1. The kagome lattice in real space (a) as described by 
its first few reciprocal lattice vectors (b) . Structure factors for 
the included reciprocal vectors are shown in black and prim- 
itive lattice vectors kp as gray arrows. The kagome lattice 
(illustrated in black, with basis and lattice vectors in white) 
is indeed determined by this limited part of the reciprocal 
lattice. 



maximum contribution to the Fourier transform p(k) = 
X^rjGbasis^*''^'^' at the reciprocal points |Gi| = ^ is 
given by G {(0, 0), (0, 1), (1/2, ^3/2), (3/2, ^3/2)}. 
With a particle density of approximately three particles 
per primitive cell any subset of three of the points vj 
is a basis for the kagome lattice. Including a hard core 
repulsion makes Vj = vi^ j ^ I, unfeasible without ener- 
getically benefiting any configuration nor affecting the 
kagome lattice. Therefore a negative perturbation of 
V{k) at /c = ^ breaks the degeneracy in favor of the 
kagome lattice which becomes the unique ground state. 
An energy spectrum constructed by following these pre- 
scriptions is shown in Fig. [2l 

A typical energy spectrum fulfilling the above con- 
straints results in a highly complex potential, not real- 
izable in experimental systems. However, the properties 
of the energy spectrum that select the ground state are 
far from a complete specification. This has dual conse- 
quences in that, on one hand, there is considerable leeway 
in varying experimental realizations for a target system 
but, on the other hand, the prescriptions leaves a fairly 
large space of possible designs to be searched. While 
having no effect on the ground states, these choices can 
greatly affect the realizability of the potential as well as 
the extent to which a system of particles interacting with 
the resulting potential can achieve the ground state un- 
der realistic dynamics. Reasonable requirements on the 
designed potentials include smoothness and restrictions 
on their range. To understand how this can be achieved 
we note that the Fourier relation in Eq. ([2]) implies that 
there exist a version of Heisenberg's uncertainty princi- 
ple relating the smoothness of the potential in real space 
to a screening in reciprocal space, and vice versa. Sev- 
eral general principles on how the power spectrum should 
be designed follow from this observation. Restricting 
the power spectrum to the low k region gives smoother 
potentials, evident in Fig. [2l Notice how the first few 
minima of the power spectrum is sufficient to determine 
the lattice periodicity. A screening of the power spec- 
trum provides further smoothing of the potential [19|. 
Conversely, to get a short ranged potential the spectrum 
should be smooth; in Fig.[2]this is achieved using third or- 
ders splines as interpolations between the extreme points. 

To test our method we perform Monte Carlo simula- 
tions with the constructed interaction potentials to con- 
firm that they reliably and in a positive density range 
self-assemble into their respective target lattices. All sim- 
ulations starts in a random disordered state and we con- 
sider only local moves of the particles in the Monte Carlo 
updates. We do not fine-tune the dimensions of the simu- 
lated systems, therefore some deformations of the lattices 
are necessary for it to fit. The simulations confirm that 
the designed potentials are robust to such complications. 

Fig. [2^ illustrates the relaxation process for the physi- 
cally most interesting structure, the kagome lattice. The 
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FIG. 2. A spectral method for targeted self-assembly. By designing an energy spectrum V{k) (b) with minima at the peaks of 
p and taking its Fourier transform we arrive at a potential V{r) (c) that causes particles to self-assemble into a target lattice, 
(a) Time evolution of a self- assembling kagome lattice. The number of grains diminishes by random walk of their boundaries, 
a process faster for smaller systems (d) but where some local defects are still present as the dimensions of the system is not 
fine tuned to fit the lattice. 



snapshots of the time evolution show how the lattice as- 
semble locally into grains with boundaries that diffuse 
and slowly annihilate. The system size in Fig. [2^ is 
chosen large enough to illustrate this process but makes 
complete relaxation too computationally demanding. In 
Fig. [2]i we show the global ground state for a smaller sys- 
tem with the same potential. It should be noted that the 
kagome lattice is significantly more complex than lattices 
previously achieved through targeted self-assembly with 
isotropic interactions. 

To demonstrate the range of our method, we also as- 
semble a snub square tiling with isosceles triangles [l6| , 
shown in Fig. O This structure is even more complex 
than the kagome lattice since it has a square Bravais lat- 
tice and requires more perturbations to single out the 
correct basis. We use perturbations at four points in 
the energy spectrum, three positive perturbations where 
p{Gi) = and one negative where p{Gi) = n. 

For completeness we also show how our method applies 
to the previously self- assembled honeycomb lattice [4]. It 
is significantly simpler, but physically interesting espe- 
cially in its graphene incarnation [20]. Since the power 
spectrum has fewer features {p{Gi) G {2,1}), even a rel- 
atively simple potential can exhibit the crucial properties 
assuring that the honeycomb lattice will self-assemble. 

All examples above utilizes key features of the lattice in 
reciprocal space, namely that the power spectrum reaches 
its theoretical maximal limit at reciprocal lattice sites 
\Gi\ = k for some k. If this criterion is not fulfilled, as is 
the case for several Archimedean tilings like the equilat- 




FIG. 3. Isotropic potentials that causes particles to self- 
assemble into honeycomb (a) and snub square (b) lattices. 



eral snub square, one cannot reduce the group of possi- 
ble basis configurations to a finite number (discrete po- 
sitions) and the choice of perturbations is not as straight 
forward. However, allowing for competing perturbations 
expands the group of achievable lattices. Since the most 
distinct features of the energy spectrum is those deter- 
mining the Bravais lattice, relevant for the interactions 
between particles of different unit cells, the method is 
simpler to apply to lattices with small unit cells. 

In general, proving that a certain configuration is the 




4 



ground state of a given potential is a very hard problem. 
In fact, the exact nature of the ground state is not rig- 
orously known even for simple interactions such as the 
Lennard- Jones potential [2l|. In this Letter we have de- 



scribed a direct method to design potentials for targeted 
self-assembly of lattices, a problem usually approached 
using iterative methods involving repeated relaxations 
of the system [l, Q . From our construction follows the 
somewhat counterintuitive observation that it is actually 
simpler to find a potential with a given configuration as 
a ground state than to determine the ground state (s) of 
a given potential. 

To conclude, we have demonstrated that there exist 
isotropic potentials that produce complex lattices such 
as the kagome. Previous work addressing targeted lattice 
formation have focused on the structure of the potential 
in real space, e.g. [H, whereas our approach focuses 
on the energy spectrum of the potential. We believe our 
method to be better suited for targeted self-assembly of 
lattices since their translational invariance is more natu- 
rally expressed in reciprocal space. 

The ability to design interactions that self-assemble 
into complex lattices is interesting in itself, but this re- 
mains a theoretical observation if the potentials are too 
compHcated to be feasible except in computer simula- 
tions. By analyzing how smoothing and screening affect 
the energy spectrum we can systematically simplify the 
potentials and thereby design interactions that should be 
possible to implemented in experimental systems. 
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